home *** CD-ROM | disk | FTP | other *** search
/ WINMX Assorted Textfiles / Ebooks.tar / Text - Mathematics - Numerical Mathematics and Computing (F).zip / rk4.f < prev    next >
Text File  |  2002-06-11  |  810b  |  35 lines
  1. C
  2. C PAGE 315-316: NUMERICAL MATHEMATICS AND COMPUTING, CHENEY/KINCAID, 1985
  3. C
  4. C FILE: RK4.FOR
  5. C
  6. C RUNGE-KUTTA METHOD OF ORDER 4 FOR SOLVING AN INITIAL VALUE PROBLEM (RK4,F)
  7. C
  8.       EXTERNAL F
  9.       DATA T/1.0/, X/2.0/, H/7.8125E-3/, NSTEP/72/
  10.       CALL RK4(F,T,X,H,NSTEP) 
  11.       STOP
  12.       END 
  13.   
  14.       FUNCTION F(T,X) 
  15.       F = 2.0+(X-T-1.0)**2  
  16.       RETURN
  17.       END 
  18.   
  19.       SUBROUTINE RK4(F,T,X,H,NSTEP)   
  20.       PRINT 3,T,X 
  21.       H2 = 0.5*H
  22.       START = T   
  23.       DO 2 K = 1,NSTEP      
  24.         F1 = H*F(T,X)       
  25.         F2 = H*F(T + H2,X + 0.5*F1)   
  26.         F3 = H*F(T + H2,X + 0.5*F2)   
  27.         F4 = H*F(T + H,X + F3)
  28.         X = X + (F1 + F2 + F2  + F3 + F3 + F4)/6.0
  29.         T = START + H*REAL(K) 
  30.         PRINT 3,T,X 
  31.    2  CONTINUE
  32.    3  FORMAT(5X,2E22.14)    
  33.       RETURN
  34.       END 
  35.